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We study single hole motion in LaMnOs using the classical approximation for JT lattice dis- 
tortions, a modified Lang-Firsov approximation for dynamical breathing-mode phonons, and the 
self-consistent Born approximation (verified by exact diagonalization) for hole-orbital-excitation 
scattering. We show that in the realistic parameter space for LaMnOa, quantum effects of electron- 
phonon interaction are small. The quasiparticle bandwidth W ~ 2.2J in the purely orbital t-J 
model. It is strikingly broadened to be of order t by strong static Jahn- Teller lattice distortions 
even when the polaronic band narrowing is taken into account. 
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LaMnOs is a typical parent compound of a class of 
colossal magnetoresistance materials i^i^x^xMnOa {R 
— La, Pr, Nd, Sm and A — Ca, Sr, Ba). As tempera- 
ture decreases, it undergoes a structural phase transition 
at Ts ~ 750 K from cubic to tetragonal because of Jahn- 
Teller (JT) lattice distortions of its MnOg octahedra. Be- 
low Tjv ^ 140 K, it is an A-type antiferromagnetic (AF) 
insulator. Staggered orbital ordering was recently ob- 
served in the ferromagentic planes of LaMnOa 0] , which 
can be driven by either the intra-atomic Coulomb interac- 
tion in the eg orbitals j^J^ or the cooperative JT splitting 



of the degenerate eg orbitals 



When electrons are 



removed from the Cg orbitals upon doping, charge fluctu- 
ations as well as spin, orbital defects are introduced into 
the system. LaMnOa with one hole is one of the sim- 
plest real systems containing charge, spin, orbital, and 
lattice dynamics. The clear understanding of single hole 
motion in it is an essential step to the full understanding 
of abundant dopant-induced phase transition phenomena 
observed in this class of materials |^ . 

The Coulomb interaction (~ 7 eV) in the eg orbitals 
eliminates doubly occupied sites and results in a two- 
dimensional anisotropic orbital t-J model [Q. There ex- 
ist off-diagonal transfer matrix elements, thus in principle 
a hole can hop without disturbing the staggered orbital 
background. If the orbitals were frozen, a hole would 
move freely and disperse with bandwidth 2t, as described 
by the large-C/ LDA-I-C/ band calculation However, 
the low-energy physics may be controlled by two polar- 
ization effects: orbital polarization and lattice polariza- 
tion. First, the motion of a hole will distort the orbital 
order via its diagonal transfer matrix elements; this will 
induce a strong polarized cloud of orbital excitations in 
the vicinity of the hole. The resultant quasiparticle (QP) 
is of orbital-polaron type with bandwidth cx J and a shift 
of spectral weight from the coherent to the incoherent 
part of the hole spectrum Second, when a hole is 
present, it attracts the surrounding oxygen ions equally, 
giving rise to a breathing distortion energy. In the sin- 
gle hole problem, charge fluctuations exactly locate at 
where the hole is and thus accompany the hole propa- 



gation. One thus needs to take into account dynamical 
breathing-mode phonons. Beyond a critical hole-phonon 
coupling strength, the hole is self-trapped in an "anti- JT" 
small polaron state via the polaronic band narrowing ef- 
fect 10 -1^. Therefore, the QP behavior is sensitive to 



model parameters. In this Letter, we present a system- 
atic study on single hole motion in LaMnOa in a variety 
of parameter regions of physical interest. 

The total Hamiltonian consider here is H = Ht-j + 
Hi,r + HjT |,g,|ll],0, where 

Ht-j^~ E (i?/^Ldj6 + H.c.) 

+ ^ E ^^i^i + ^'^i'^i' T %/3(Il^r/ + T-^T^)], (1) 

(ij>ll 

i/br = -y/gbr^o y^nf (gt + O;) + cjq a\a-^, 
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where d 



t -^t 



d\^{\ ~ nia) is the constrained fermion op- 

(dld.^ - 



erator for the Cg electron at orbital a. T( 



dI^d^)/2 and 7]^ = (dl^d;^ + dI^di^)/2 are orbital- 



pseudospin operators with | |) 



d^ 



and I i) 



The anisotropic transfer matrix elements are 



3t/4, til = t/4, and tl> = t^J = T\/3i/4, here the 
=p sign distinguishes hopping along the x and y directions. 
The orbital superexchange J-induced orbital ordering re- 
produces the experimental observation: (| T) + I i))/>/2 
and (I T) ^ I i))/V2 in the A and B sublattices respec- 
tively I],!!. In i/br, "i' = 1 - Z^a dladia is the hole num- 
ber operator, ai's are the breathing-mode phonon opera- 
tors with frequency loq- -E-br is the hole-phonon coupling 
strength. In Hj^ |l^, Q2i and Qsi are, respectively, JT 
distortions for the 3z^ — and x'^ — modes satisfying 
Q2i = gcosV'i and Qsi = gsin-^i with -^igA = -V'je-B 
iQ. For simplicity, V'ieA is set to be 7r/2 so that the 
orbital ordering driven by the cooperative JT effect also 
agrees with the experiment. 

We employ the slave-fermion formalism to cope with 
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the constraint of no doubly occupancy |T^. First, we 
perform a uniform rotation of^rbitals by 90° about the 

—Tf, Tl^ ^ T? in order to obtain the Neel configura- 

tion I • f+^rr+s---) = I-- - itit •••>• Then, 

considering the Neel state as the vacuum state, we de- 
fine holon (spinless fermion) operators hi so that Ci| = 



Cix = hi on the | sublattice and Cj| 



hj on the f sublattice. Here h = T ^ on the | sublattice 

and on the t sublattice are hard-core boson opera- 
tors. Such a treatment is in spirit similar to that done 
by Schmitt-Rink et al. for one hole motion in a quantum 
antifcrromagnet |r^-|l7|]. 

Following Roder, Zang, and Bishop we 
treat dynamical phonons within the modified Lang- 
Firsov approximation using a canonical transforma- 
tion H = U^HU, U = e-^if'^^e-'^^^^), where 
5i(A) — A/2^/Eh^u!o ~ ^i) with A measur- 

ing static breathing-mode distortions and S2{t) = 
—T^ E^iijujQ 'Ylii hlh^{al — Cj) with r measuring the de- 
gree of the polaron effect. Averaging H over the phonon 
vacuum, we arrive at an effective orbital-lattice-polaron 
Hamiltonian H = Eq{t) + Hcs 

k q 

+ E ^ k/lk-q(Mkq/3q + ^kq/^q+q) + H.C. (2) 
kq 

where Q = (7r,7r). Eo{t) = Nx[il - r)A + A^/4,Eb, - 
£'br(2T — T^)] with X ~ 1/N being the hole concen- 
tration and A — 2£'br(l — t)x. The hole-orbital- wave 
coupfing functions are Mkq = ^(7kUq + 7k-qWq), 

^kq = -^(»ykWq-'7k-q'"q) with ?iq = {[{Aq + Bq) / UJq + 
l]/2}l/2 and «q = -Sgn(i?q){[(ylq + i?q)/c^q - l]/2}l/2. 

Here short notations are 7k — (coskx + cosky)/2, rj^ = 
{cosk^ — cos ky)/ 2, Aq = 3J+i?jT, and =J7q/2. The 



with 



/Sq's are orbital wave operators, bq 

dispersion ujq = Aq{Aq + 2Bq). The bare hole disper- 
sion is S]i{T) — — C(T)t7k with ^(r 



exp(-£;brT^/wo) 



being the polaronic band narrowing, ^(t) is determined 
by the following equation [n2|: 



Eu 



-ujQ In^ 



1 



dEn^iniO 



(1 + x)u;o 



(3) 



where i?inin(C) is the minimum of the QP dispersion i?k = 
Ek-f Rel](k, Ek) with S(k, oj) being the holon self-energy. 
Note that in this effective Hamiltonian, the JT effect acts 
as EjT = 2\^/K in Aq. 

We calculate the holon Green's function G'(k, w) = [lo— 
Ek — S(k, cj)]^^ treating the hole-orbital-wave coupling 
within the self-consistent Born approximation (SCBA) 
. The self-energy is thus of the form 
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FIG. 1. QP spectral functions A(k,ui) — — ImG(k, a;)/7r for 
the orbital t-J model on the 4x4 cluster calculated using the 
SCBA (left panel) and using the ED technique (right panel). 
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FIG. 2. QP dispersion for the orbital t-J model with 
J — 0.3t: The ED results for the full Hamiltonian (solid 
squares) and the truncated Hamiltonian (solid circles); The 
SCBA results on the 4x4 cluster (open circles) and the 
20 X 20 one (open squares). Solid lines are the fits using 
eI = ao + ai7k + 02 cos cos ky + a372k. 



S(k, U;)^Y1 [^kqG(k - q, - UJq) 

q 

+ ^kqG(k-q,C^-Wq+Q)] 



(4) 



We test the applicability of the SCBA by numerically 
diagonalizing Ht-j on 16- and 18-site clusters. In ad- 
dition, since in the SCBA the AF orbital background 
is treated within linear orbital-wave theory, the mixed 
terms cx T^T? + T.^T? in i?t-j, which contribute only 
in higher order orbital-wave theory, are neglected in the 
SCBA. It is thus interesting to numerically diagonalize 
a "truncated" Hamiltonian without these terms ||^. As 
shown in Fig. 1 and Fig. 2, we find that all of the SCBA 
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results are in good agreement with the exact diagonahza- 
tion (ED) results, especially with those for the truncated 
Hamiltonian as expected. 

The physically relevant parameter space for LaMnOs 
are: t ~ 0.72 eV (energy unit) and 0.1 < J < 0.3 are es- 
timated from photoemission experiments ||l8|. ujq = 0.1 
is taken from Raman data |jl^. Ejt = 1.5 and i^br is 
of order Ejt 111^. Numerical studies found that 
the A-type spin and C-type orbital structures were sta- 
bilized in this region of parameter space Q . Below let us 
discuss the QP properties in different parameter regions 
of interest. All calculations are performed on the 16 x 16 
lattice unless noted. 

First, we consider the single hole motion in the pure 
orbital t-J model (i.e. i?br = -Ejt = 0). As displayed in 
Fig. 1, at any k, there is a well-defined quasiparticle pole 
(i.e., zero orbital- wave) at the low energy side which is 
well separated from a broad, incoherent, multiple-orbital- 
wave background extending to the full free-electron band- 
width. The bottom and the top of the quasiparticle (QP) 
band locate at (0,0) and (tt, tt), respectively. In Fig. 
2, the spectrum of the orbital polaron is flat at large 
momenta, which leads to a strongly distorted density of 
states with a massive peak at the top of the QP band. 
In the realistic range of 0.01 < J < 0.5, the QP band- 
width W ~ 2.2J scales with J (see Fig. 3). This new 
low energy scale is quite similar to that in the cuprate t-J 
model where the spin is conserved during hole hopping 
|p5|-p7t. In the latter, the staggered spin background 
is disturbed by hole propagation, leading to a vanish- 
ing bare hole dispersion e^, and thereafter restored by 
Heisenberg spin flipping, forming a QP band with width 
~ 2 J (the characteristic energy of spin waves). In the 
present case of Ek — —tj-k, the energy scale of ~ 2.2 J 
can be understood in the following way: For small J, 
orbital excitations are easily stimulated by incoherent 
hole motion and accompany the hole propagation. Thus 
the hole spectral weights Z(k) = [1 — d'E(k,uj)/duj]^^^^ 
are strongly reduced by such a cloud of polarized orbital 
waves. The reduced weights by incoherent hole motion 
can be approximately obtained by neglecting in the 
calculation, referred to as ZgCk). Then one can naively 
express the QP dispersion as Zo(k)ek- In the range of 
0.01 < J < 0.5, ^o(k) ~ l.3J/t. Therefore, W scales 
with J. 

Second, we consider the case of strong static JT distor- 
tions. The contribution of static JT interaction is adding 
an Ising-like component to the excitations and inducing 
a large gap in the orbital excitation spectrum. Thus, the 
JT effect stabilizes the orbital ordering. As a result, even 
for small J, orbital excitations are difhcultly stimulated. 
In fact, for ^ £,{T)t due to either large i?jT or large 
i^bri Eq. (^ can be solved analytically in perturbation 
theory: ~ -C(T)t7k - 0{^H^/Aq). Thus, the QP 
band is narrowed mainly by latttice polarons instead of 
orbital polarons. The QP spectral weight is 



W/2t 




FIG. 3. QP bandwidth W as a function of J at iSbr ~ 0. 




FIG. 4. The variation of bandwidth W and the polaronic 
narrowing ^(r) with Ehr at J = 0.3 for (a) lj = 0.1 and (b) 
uj = 0.5. 

^{t) - 0(C^i^M2) approximately. The QP bandwidth 
scales with t. As illustrated in Fig. 3, > 1.5t is a 
slowly changing function of J at i?jx ~ 1.5. Fig. 4(a) 
shows W and ^(r) as a function of i?br at i?jT = 1.5 and 
luq = 0-1; they behave similarly in general, indicating 
that lattice polarization overwhelms orbital polarization 
in a wide range of i^br- For not too large i^br, the small 
deviation of W/2t from ^(r) is attributed to the remnant 
hole-orbital-wave coupling. W decreases exponentially 
with Ehr for Ehr > 12.5. This is the Lang-Firsov band 
narrowing. However, for the realistic value of i^br < 5, 
W > 1.2t. Therefore, in the presence of strong static JT 
distortions, the QP bandwidth is strikingly broadened in 
comparison with W ~ 2.2 J obtained in the purely orbital 
t-J model even when the quantum effect of polaronic 
band narrowing is taken into account. 

Third, it is interesting to examine the case of no static 
JT distortions (i.e., i?jT = 0), which is relevant to the 
experimental fact that the static JT distortion rapidly 
disappears around a; ~ 0.1 in Lai_j;Srj;Mn03 [Q. For 
not too large i?brj there is mix of polarization of phonons 
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and orbital-waves. Fig. 4(a) shows the QP bandwidth 
W as a function of i?br at ojq = 0.1. For Ejt = 0, 
W/2t remains almost unchanged I.IJ) in the range of 
< i^br < 10 and obviously deviates from the polaronic 
band narrowing ^(t) up to i?br — 17. This implies that 
the QP is of orbital-polaron type for _Ebr < 10 and is of 
polaron type for Sbr > 17 as well as a mix of the two 
types for 10 < E\,^ < 17. Therefore, for realistic value of 
Sbr < 5, the electron-phonon interaction effect on the QP 
band can be neglected. The hole-orbital-wave scattering 
dominates the formation of the QP. 

Finally, it should be made clear that the small value of 
ciJo 0.1 (which is however the highest value of phonon 
frequency observed in Raman spectra |l9[| ) underlies the 
unimportance of quantum effects of electron-phonon in- 
teraction on the QP band. Note that the lattice polaron 
effect is the most pronounced in the antiadiabatic limit 
(luq — + oo), where 7 — > 1, and is negligible in the adia- 
batic limit {ujq 0), where 7^0. At o^o = 0.1, r < 0.1, 
the system is close to adiabatic limit. Fig. 4(b) shows 
the same quantities as Fig 4(a) but at ujq = 0.5 which is 
chosen in Ref. |12|. The quantum effects of hole- phonon 
interaction are improved remarkably. For example, the 
critical values of E'br at which different QP types occur 
are 3 times as smaller as those at loq = 0.1. 

Summarizing, we present a systematic study on single 
hole motion in LaMnOs. We show that in the realistic 
parameter space for LaMnOa, quantum effects of hole- 
phonon interaction are small. In the purely orbital t-J 
model, the QP bandwidth W scales with J. Considering 
the hole-phonon coupling, W, which however scales with 
t, is strikingly broadened by strong static JT lattice dis- 
tortions even when the polaronic band narrowing is taken 
into account. We predict that orbital polarization is pro- 
nounced for X > 0.1 where static JT distortions disap- 
pear. Our results can be tested by future angle-resolved 
photoemission spectroscopy (ARPES) experiments. 

The authors thank P. W. Leung, F. C. Zhang, and H. 
Zheng for useful discussions. The ED program was imple- 
mented on the basis of C codes provided by P. W. Leung. 
W.G.Y. is supported by CUHK 4288/OOP 2160148. 

Note added. — Very recently, two papers concerning 
ARPES of LaMnOa have appeared. Perebeinos and 
Allen addressed the Frank-Condon broadening ef- 
fect driven by electron-phonon interaction on the bare 
hole dispersion in the JT-ordered ground state. Brink, 
Horsch, and Oles (BHD) calculated the single- hole 
spectral functions for the purely orbital t-J model using 
the SCBA and discussed the crystal-field effect on the QP 
band. Their pieces of work are in a sense complementary 
to ours. The present work takes into account both of or- 
bital polarization and lattice polarization, giving rise to 
a different, yet more realistic picture from BHO's. 
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